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Spatial pattern formation in excitable fluctuating media was researched analyti- 
cally from the point of view of the order parameters concept. The reaction-diffusion 
system in external noise is considered as a model of such medium. Stochastic equa- 
tions for the unstable mode amplitudes (order parameters), dispersion equations for 
the unstable mode averaged amplitudes, and the Fokker-Planck equation for the 
order parameters have been obtained. The developed theory makes it possible to 
analyze different noise-induced effects, including the variation of boundaries of or- 
dering and disordering phase transitions depending on the parameters of external 
noise 
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I. INTRODUCTION 

Noise is present in real systems of any type. The influence of external noise on nonlinear 
open distributed systems is very diverse and sometimes leads to quite unexpected results. It 
is known, that noise-induced spatial patterns and fronts [§1, @], noise-induced resonant 
pattern and frequency- locking phenomena [101 ] . pure noise-induced phase transitions lllk 

noise-induced spatiotemporal intermittency |l6| . 
noise-supported traveling and noise-sustained 



141 ]. noise- induced phase separation 



spatiotemporal stochastic resonance 



15 



17-19 



convective structures [2~ol-22]. noise-induced synchronization (23-25], etc., may arise in such 
systems. 

Theoretical study of spatio-temporal dynamics of nonlinear open distributed systems is 
carried out using different methods 

Stability of the homogeneous state with respect to small perturbations can be analyzed 
in linear approximation. However, linear approximation is not acceptable to describe the 
evolution of the system near the threshold of self-organization, as the phenomenon self- 
organization itself is essentially a nonlinear effect. 

There is an approach to the study of noise-induced phenomena, based on the well-known 



mean- field approximation |4|, |12|, |15|, |26|, |27| . In this approximation it is assumed, that the 
interaction between a certain spatial point and its nearest neighbors occurs through the field, 
whose value corresponds to the statistically average field at this point. Herewith, a suitable 
way is used to carry out the discretization of the space of the initial distributed system 
and the Fokker-Planck equation (FPE) for the multivariate probability density function can 



2 



be written for field values in the points received by a regular lattice. The obtained FPE 
is integrated over the values of the field at all points except the given one. This leads to 
FPE for the one- dimensional probability density function values of the field at a given point. 
In the latter equation the conditional average values of the field at neighboring points are 
replaced by an average value of the field at a given point. This approach can predict the 
existence of noise-induced "disorder-order-disorder" phase transitions. 

Another analytical approach to the study of nonequilibrium phase transitions with the 
spatial pattern formation is based on generalized Ginzburg-Landau equations for the order 
parameters of the systems. This approach makes it possible to describe the behavior of the 
system near the transition point. It is based on the separation of eigenmodes of the system 
into damped (stable) and undamped (unstable) modes (order parameters), and into adia- 
batic elimination of stable modes. The method of generalized Ginzburg-Landau equations 
for the systems, the right-hand side of which contains additive white noise, is described in 

In papers [14|, |34[ spatio-temporal evolution of the nonequilibrium extended systems is 
investigated by means of dynamic renormalization groups. In [l4| it is shown that under 
certain conditions a new genuine nonequilibrium universality class arises due to the presence 
of multiplicative noise. 

Moreover, there are approaches based on anal ysis of correlation functions of the dynamic 
variables of the system or structure functions j3l.]l~2"j]. on the study of higher-order moments 
' 13, and others 0, M, 0, Q . 



The aim of this paper is to develop a theory, which would allow us from the unified point 
of view of the concept of order parameters, to carry out a consistent and detailed study of 
spatial pattern formation spontaneously arising in open nonlinear distributed systems with 
external noise both in the vicinity of the transition point and away from it. 



II. STOCHASTIC EQUATIONS FOR THE ORDER PARAMETERS 

Systems of the reaction-diffusion type 

— = P fc (xi, x 2 , x 3 , 77, r, t) + D k V 2 x k , k = 1, 2, 3 . . . (1) 

are one of the basic models of a nonequilibrium excitable medium. In equation (pQ) x k are the 
medium state functions (dynamic variables), Pk(xi, x 2 , £3, 77, r, t) are nonlinear functions, 
that define the interaction and the evolution of the component x k in space and in time, 
Dk are diffusion coefficients of components, 77 = (%i,. . . ,Xn,Vi,- ■ ■ Ms) is the vector, whose 
components are the control parameters describing the effect of the external environment 
on the system. Further, without loss of generality, we consider a two-component system of 
type (0Q). Nevertheless, the proposed research method is easily extended to multicomponent 
systems. 

In external noise the spatio-temporal dynamics of nonequilibrium systems ([1]) for k = 2 
can be described by the following system of equations: 

a m 

-J£L = p 1 ( Xl ,x 2 ,Xio, -,XmO, -,Xn) + E flj (r, t)P\j (Xi , X 2 , Xm+1 , -,Xn') + £>lV 2 Xi, 

iff = i^Ei, 2:2,7710, ...,7/10, :.,r]s) + E f 2 j(r,t)P 2 j(x 1 , x 2 ,r]i +1 , ...,rf a >) + D 2 W 2 x 2} 

i=i 
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where m and / are the number of fluctuating parameters in the first and second equations, 
respectively, Xjo, Vjo are spatio-temporal average parameters, fij(r,t) (i=l,2) are random 
fields describing the noise of the appropriate parameters with respect to their mean values 
with(/ o -(r,t)) = 0. 

We define the statistical properties of random fields fij(r,t) according to the properties 
of the environment. Fluctuations in the environment represent the summarized effect of 
many weakly coupled factors. It follows from the central limit theorem, that fluctuations 
of the external source have a Gaussian distribution. The ergodic Markovian and Gaussian 
properties of fluctuating environment limits the choice of random fields for modeling the 
fluctuations of the environment by a stationary homogeneous isotropic Gaussian field with 



the exponential time- correlation function [28 



#[ii i (r,^),^ i^ (r / ,t , )] = $ l (|r-r , |)exp(-A; t4 |t-t / |)^^v, (3) 

where $j(|r — r'|) define the spatial dependence of correlation functions of the random fields. 
The cross-correlation of the fields /^(r, t) and fiij/(r,t) is absent. The correlation time 
r t = k^ 1 is the characteristic time scale of external noise. Hereafter, we use the notation 
K[F 1 ,F 2 ], that is defined by the equality K[F 1 ,F 2 ) = (FiF 2 ) - (F 1 )(F 2 ) for the correlation 
function. 

Let us assume for simplicity, that m — I — 1, and introduce the dimensionless variables 
T — Xwt an d r' = r^Jxio/ F>i, where Xw is a parameter, that has the dimension of inverse 
time. Now the system (J2]) can be rewritten as 

= P{(x 1 ,x 2 ,Xio,-,Xn) + fn(r',T)P[ 1 (xi,x 2 ,X2,--,Xn>) + V ,2 xi, ^ 
§^ = P2(xuX2,Viq,-,Vb) + hi(r',T)P^(x h x 2 ,rj 2 , ...,rj gl ) + DV 2 x 2 . 

Functions with primes are different from the corresponding functions without primes by a 
multiplier xlo- D = D 2 /D 1 . Hereafter, the primes are dropped for simplicity. 

Suppose, that in the deterministic case the values of parameter Xii- ■ ■ iXniVii- ■ ■ jVs 
are such, that there are stable stationary states xio and x 2 o defined by the equations 
Pi(xi, x 2 , xi, Xn) = and P 2 (xi, x 2 , r?i, ...,T] S ) = 0. 

We write the equations (jl]) in the operator form. Simultaneously, we select from its 
right-hand side linear i^(V 2 )q, nonlinear g(q), and random F(r ; r) components |29| 

^-K(V 2 )q = g(q) + F(r,r). (5) 

Vector q describes the deviation of the dynamic variables from their equilibrium values: 
q = (xi — xio,x 2 — x 2 o). A linear operator K(V 2 ) takes the form 

, 2 . ( a n + V 2 a 12 \ dp 



a 21 a 22 + DV 1 ) ' 11 dxj 



, i,j = 1,2. (6) 

Vector g(q) contains quadratic and cubic nonlinearities obtained by series expansion of 
the right-hand deterministic side of the equation (jl]). Its components are defined as follows 

2 2 

&(q) = 9i!luQ^Qu + otluA^q^ (7) 



where a {2) = ± d2p > 



fJ,,U=l fJ,,U,K=l 

(3) _ 1 d 3 P t I 
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Vector F contains the random fields: 

F=( f 11 ( T > T ) Pi - 1 ( Xl > x *>X 2 >"->X«>')\ 

\f2l(r,T)P 2 l(x 1 ,X 2 ,T]2,-,Vs>) )' 

To research the stability of the stationary state of a deterministic system, we assume 
that the vector q has the form q = q exp(Ar + ikr). The respective characteristic equation 
A 2 — aX + = has solutions 



A l,2l k J - , (SJ 

where a(k) = Trfa) - (1 + D)k 2 , (3(k) = Debfa,) - (Da u + a 22 )k 2 + Dk 4 . 

Conditions a(k) < and /3(k) < define aperiodic instability, herewith one positive real 
root of eq. (EJ) appears: Re(X\) > 0, Im(Xi t2 ) = 0. Further we shall consider only this case. 

Suppose, that one of the parameters xiv • • ,Xn,Vh- ■ ■ ils, for example %2, is bifurcation, 
i.e., at some critical value of this parameter \2c there exists a critical wave number k c , when 
the conditions Re{Xi{k c )) = 0, = are fulfilled. 



dk 



k=kc 



Represent the vector q in the form of superposition 

q(r,r) = EO (j) (k')#(r)e jk ' r , (9) 

where O^(k) are eigenvectors of operator K(V 2 ), $ (r) are unknown amplitudes, exp(zkr) 
are eigenfunctions of operator V 2 . Here it is assumed, that the vector q(r,r) is a superpo- 
sition of plane waves, but depending on the symmetry of the problem Bessel functions or 
spherical wave functions are to be chosen as the eigenfunctions. 

The vector F contains nonlinear functions Pa. They need to be expanded in a Taylor 
series in powers of the components of vector q. We restrict our consideration to quadratic 
terms in the expansion of Pa. It is easily shown, that the terms of this order are necessary 
to obtain the dispersion equation for the averaged amplitudes of the unstable modes with 
an accuracy to terms quadratic in the fluctuation intensity. As a result, the components of 
vector F containing random fields take the form: 

2 2 

Fi = /ii(r,r) [ Pi 0) + ^2pQq fl + J2 pf^nQv 

fl=l fl,U = l 



J2) = 1 8*Pu 
2! dxndx v 



where pf> = Pii{x w ,x 20 ), p$ = 

Unstable modes lie in a narrow band of wavenumbers defining the instability region of 
the system. This makes it possible to construct wave packets by summing over the wave 
vectors, which lie in a small interval. Thus, carrying modes with discrete values of wave 
vectors and slowly varying amplitudes £k^( T ) are chosen [29J. 

To obtain the equations for the mode amplitudes £jf (r) we substitute eq. (jSJ) in eq. (j5J), 
multiply the equation derived on the left by exp(— zkr)0*^ ^ (k) and integrate over the region, 
which is much greater than the oscillation period exp(ikr), but in which ( T ) varies very 
little. Here 0*^'^(k) are eigenvectors of the operator conjugate to K(V 2 ): OV'O*^) = 

After transformations, the system of equations for the amplitudes of the modes takes the 
form 
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%-A,(k)^ } = e ^,,(k ) k^k'')^;' ) ^: ) ^(k / +k''-k)+ (io) 

j'i",k'k" 

+ E ^"i"'( k > k '' k "' k " , )^" ) e&" ) e k C ) 5(k' + k" + k'" - k) + E 0? ) (k)pW^(r)+ 



i'i"i"',k'k"k"' ip=l 



+ EE <M k > k ')^ %,k-*(r) + E E 4^"( k > k', k")# } & %, k - k - k »(r) 

cp=lj',W <p=lj'j",W 

Functions o$„(k, k', k"), ^^(k, k', k", k'"), ^( k ' k ')' £ S'i"( k > k '> k ") introduced in 
equations (TTOT) are presented in the Appendix A. 

Random processes z^^ir) = J /^i(r, r)e~ lkr <ir are components of the random vector field 
z(t) with zero mean, tp and k are index arguments of this field. 

Assuming, that the correlation time of random fields is considerably smaller than all the 
characteristic times of the deterministic problem fl4]), correlation functions for the compo- 
nents of field z(t) will have the form: K[zj^(t), zi^(t)] — gji(\k\)S(k— k')S(t— r)5ji. Further 
we assume for definiteness 

$ J (|r-r'|) = ^exp(-A; / ,|r-r'|). (11) 

Here 9j are noise intensity, kfj are magnitudes inverse to correlation length. For two- 
dimensional media g^ = 9jkfj/[2n 2 (k 2 + ~ 3 / 2 ]. 

The system f JTOj) contains both stable and unstable modes. In the vicinity of the bi- 
furcation point the relaxation times of unstable modes are considerably greater then the 
relaxation times of stable modes, therefore the latter adiabatically follow the former ones. 
This makes it possible to exclude the stable modes from the equations f lTU]) . To perform 



the procedure of adiabatic elimination [29|, |30| of stable modes, we rewrite the system of 



equations (fTOj) dividing it into two subsystems of equations for the unstable modes (denote 
them by an additional index (u)) and for the stable modes (s). 

Since the unstable modes can grow to infinity if we neglect the nonlinear terms, we write 
equations for them with an accuracy of cubic terms, which provide nonlinear stabilization 
of the instability (if the cubic terms are not enough, it is necessary to take into account the 
fifth orders). 

We assume that the amplitudes of stable modes are significantly smaller than the ampli- 
tudes of unstable modes |£ s | << and their variations occur self-consistently: £ s ~ In 
addition, ~ £ s . 

In the equations for the stable mode amplitudes we keep only the terms necessary to 
obtain the equations for the unstable mode amplitudes with an accuracy to third-order 
terms. Then from eq. ([TO]) for the stable mode amplitudes we obtain the following equation: 



^ " A.ikJC' = E 4 } (k s , k' u , k" u )£l$J(k> u + k" u - k s ) + (12) 



dr 

k'«k"?i 



ip=\ i/?=l k'w 



+ E 0^\k s )pfz^ s {r) + E E ^i(k s , k , u )e£^,ks-k'«(r) + 



2 

+ E E £ </,ll(ks> k u , k tj£ k / u £ k // u £ (/ , l ks-k'M-k''u(' r )- 
<fi=l k'?ik"« 



The equations for the unstable mode amplitudes have the form: 



dr 



k'«k"u 

4- V rr (1) rk k' k" k'" U^M 1 ) z^ 1 ) A/V _u k" -4- k'" - k ) + 
k'uk"wk"'u 



£ k ^ from them and substituting the latter in eq. ffT3l) . we obtain a system of equations for 
the unstable mode amplitudes £ k ^: 



V=l 1/3=1 k'u 

2 

+ E E 

(r)+ 

V=l k'uk"u 

+ E E^i(k u ,k s )e k t ) ^k u -k s (r) + 

■ip,ip=l ks 

2 

+ E E (k«) k uj k s ) + cr^ 1 (k u , k s , k u ) £ k / u £ k t <Kk « + k s — k u ) + 

V>=1 k'uks 
2 

+ E E [ £ ¥>,i^(kuj k u , k s ) + e^ 1 (k u , k s , k u ) £ k 'i£kt ^kw-k'u-ksCr)- 

i/j,</j=l k'uks 

Neglecting the time derivative d^} /dr in equations ( |T2l) [23], expressing the amplitudes 



dr 



(14) 



F to (r) = AiW^ + E ^| } (k u ,k' u ,k\)^k^(k , u + k'' ti -k u )+ 

k'jik"u 

+ E/ w (k u , k u , k u , k u , k u — k u)£k'«£k"'u£k'"'!i^(k « + k u + k u — k u ) + 

k'Mk"uk"'u 

2 2 

+Eo; (1) (k u ) P (°)z,, ku (r)- x: EW(k kg) ^ku-kg ( r ) z ip' ,ks ( r ) + 

^=1 ip,tp,ip'=l ks 



2 2 

+ E E %>(ku> k u )^, u 2; V p,kti-k'u(T) + E E ^(ku, k u , k U )^ k / U ^ k // M 2 ¥ , ]ku _ k / U _ k » u (r) 

</p=l k'u <p=l k'uk"ti 

2 

^ ] ^ ] -^■tp,ip,ifi i (k u , k s , k u^^]^ u Zipt fcs—\L l u(j'^Ztp,Vu—\is(j~) 

ip,tp,ip'=l ks,k'« 
2 

— X] ^tp,ip,ip' (k u , k s , k u^k'v^ku— k'u.— kg^^^ks^) 

ip,ip,ip'=l k'uks 

2 

— E/ E/ C^,i/>,¥/(k u , k s , k u , k ^£ k / u £ k // u £< /3 ' l ks-k'u-k''u(' r )^i 1 o,ku-ks(' r ) — 

ip,<p,ip'=l k'uk"uks 
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2 

^ l 4>,i>,f' (k-uj k u , k s , k u)^i c 'M^k"'u^,ku-k'M-ks(' r )^i ( 3',ks-k"M(' ir ) — 

ip,ip' ,ip=l k'wk"itks 

2 

— Etf )( p{k u , k u , — k u , k u , k u )^i c 'M l Ck"M'Ck'"«^,kM-k'u-k"M-k'"«(' r ) — 

i/},ip=l k'?ik"Mk"'« 

2 

— 51 ^V,^(ku, k u , k u , k u)£ k 'u£k''?i£k'''u%>,k«-k'«-k''u-k' / '<ii('7~) — 

tp,ip=l k'Mk"uk"'« 

2 111 

Gip,tp,tp'(ky,i k u , k s , k u, k u)^] t /. lt ^ k // u ^j c /// u ^( 1 D,kM-k'u-ks(' r )^(p',ks-k"u-k'"M(' r ) • 

ip,(p',ip=l k'uk"Mk"'Mks 

Function u(k u , k' u , k" u , k'" u , k u — k' u ) and others introduced in equations ( fT4l) are presented 
in the Appendix B. The modes are order parameters. Their collaboration or competition 
determines the behavior of the system. 

The system of equations (TH1) is still difficult to analyze because it contains random 
components. Further analysis of the equations (jl4jl may consist in their statistical averaging 
or in obtaining the Fokker-Planck equation. 



III. STATISTICAL AVERAGING 

For statistical averaging we used the relationship between the moments and the cor- 



relation functions [37] and mult i- dimensional generalization of Furutsu - Novikov formula 
3jJ . Taking into account the formal solution of Eqs. ( JT41) and assuming, that we can ne- 
glect higher than second order correlation functions, it can be shown, that in Eqs. (fl4l) the 

terms containing the product d'kk»ld'''^,k 1 (^)v,k 2 ( r ) and d'ld'l£k»'^,k( r ) must be 
discarded since they do not contribute to the average values of the modes £ k ^ in averaging. 
We note here, that the procedure of correlation splitting leads to the appearance of a similar 
correlation for the other interacting modes. Therefore, this procedure should be performed 
until all the terms containing the necessary degree of intensity fluctuations are taken into 
account. The remaining correlation functions can be neglected owing to their smallness, as 
the terms obtained after their splitting are proportional to a higher degree of noise intensity. 

In order to obtain after averaging corrections to the increments of unstable mode am- 
plitudes with an accuracy to terms quadratic in the noise intensity, when calculating the 
functional derivatives in the Furutsu - Novikov formula it is necessary to retain the terms 
containing the product ^ k !i^ >k (r). 

Below is given only the structure of equations obtained by the procedure of averaging 
since they have a very complicated form: 

- Ai(k u )(^) = L (K) + Li (K, k' u ) + (15) 

+L 2 (k u , k' u , k" u , k"' u , k s ) (£k''u£k'"u) + L3 (k u , k' u , k" u , k"' u , k^ y , k s , k' s ^ (Ck''wCk"'«Ck^ u )- 

The analysis of the equations ( ITS]) leads to the following conclusions. 

First, after the averaging of Eqs. (THI) additional terms, that do not depend on ^ k ^, 
arise in ffT5l) . They are determined by the parameters of the problem, the type of the 
correlation function gjj, the noise intensity, and the wavenumber of a given mode. 
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Second, in the system f|T5l) there are additional terms proportional to This leads to 

a variation of eigenvalues of unstable mode amplitudes in comparison with the deterministic 
case. As a result, the region of instability of the system, the beginning of the process of 
destruction of a statistically stationary homogeneous state and pattern formation as well as 
the duration of the transitional regime from one statistically stationary state to another are 
changed. 

Selecting from Li (k u , k' u ) the terms that contribute to the increment of (£ku\ we obtain 
the dispersion equation 

A = Ax(k u ) + L x (k„, k u ) + L a (k u , -k u ) . (16) 

When deriving eq.f fT6|) we take into account the, fact that £1^ = £^1 = £ku> because the 
solutions of equations (TH|) must be real. From equation f fl6|) and expressions Li (k n , k' u ) 
and gu, it follows, that increments (Re A) of the unstable mode averaged amplitudes are 
proportional to the intensity of noise and depend on the correlation length. Herewith, the 
intensity of noise becomes another bifurcation parameter and patterns begin to form, when 
the parameter \ie is different from \2c- Research of particular systems shows, that the value 
of the bifurcation parameter \2 is shifted to the subcritical region for the case of deterministic 
description. Finally, we note, that in the external random fields pattern formation occurs 
due to multimode interactions, whereby the conditions the resonant interaction of the modes 
are also different from the deterministic case. In particular, interaction between different 
configurations of modes with the modes with doubled wavenumbers arises, for example 
2k' n — k" u = k n , as well as the five-mode interaction. 

Thus, the above theoretical analysis allows us to describe the evolution of the stochastic 
systems under consideration, in the vicinity of the Turing bifurcation point in more detail. 



IV. FOKKER-PLANCK EQUATION FOR THE ORDER PARAMETERS 

If the system parameters are such that it is in the supercritical region then with increasing 
noise intensity, the system will "go" farther and farther from the bifurcation point of a 
deterministic system. To describe the state of the system in this case we can use the Fokker- 
Planck equation. 

For the system ( I14p the Fokker - Planck equation can be written in general form as follows 
[35]: 



F^ u {t')}dt' \w + (17) 



|^ j (J K[F ku (r),F qu (t')K'j uj. 



Here w ({£k«}' r ) ^ s multivariate probability density, which determines the probability of 
some configuration of unstable modes |£kL|- After transformation with an accuracy to 



terms linear in the noise intensity we can obtain the correlation functions appearing in (117]) : 
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F v (t)] = E Vv (K, cLu)0^\q u )p^g w (\K - q u |)5 ku _ qU)qu 5(r - 0+ 

+E M k 

a] Qui k u 2q u ) + z/(p(k u , k u 

+^77 ¥ ,(k u ,q u )77 (/ ,(q u ,2q u - kj^flku - q u \)^2du-kJ(T - t') + 

v 

+ E %>( k «> q u )^(q u , q'„, 2q u - k« - q'J# w (|k 

« q« 

<p,q'u 

+ E M k «, q«, k u - 2q u + q' u + q" u ) + z^(k u , k u - 2q u + q' u + q" tt , q tt )] x 

<^3,q'u,q"u 

xf v (q«, q' u , q /, „)^ w (|q« - q'„ - q ,, J)&Mq''«£ku-2qu+q'«+q"u£( r - 
+ 2 ^(q 

2q« + q'u, q«) + v<p(K, q«, K - 2q u + q'J] x 

tp, q'u 

X 3^(|qtl - q'J)£q'«6c? 1 -2q«+q'u ( K r ~ > 

^[F ku (r),F qu (t')] = E [0; (1) (k„)] 2 (^ 0) )% w (|k.|)5 kw 5(r - t') + 

+J2v<p{K,K - q u )O; (1) (q u )^ ^ w (|q u |)e ku -q U 5(r - 0+ 
+ E *V( k u, k '«' k « _ k '« - q M )O* (1) (q u )^ 0) 5' w (|q M |)^k'^k«-k'«-q«5(i" - 0+ 

ip,k'u 

+E r 7v( c l«> c l« ~~ k «)°? 1) ( k «H° ) 5 , w(l k «l)^q«-k«5(i" - 0+ 
+ E ^( k u, k u - q u + q' u )r/ (/ ,(q u , q'J# w (|q„ - q'J^q'u&u-qu+q'u^T - t') + 

ip, q'u 

^ ] ^(k-uj k u , k u k u q, u + q u 

)^(q«,q'«)5 , W (|q«-q / U |)Ck'«Cq'^k«-k'a-q«+q'a5(r-t / ) + 

ij5,k'u,q'u 

+ E q'u, q« - q'u - k u )O* (1) (k u )pJ 3 0) c/ w (|k M |)^ q /. [1 ^ qu _ q / u _ ku 5(r - t')+ 



+ E ^(q«> q'u> q« ~~ q'u - k u + k 'u)^(k„, k' u )5( w (|k u -k / ?1 |)^ k / u ^ q / u ^ qu _ q / u _ ku+k / u (5(r-t / ) + 

i^k'Ujq'u 

+ E ^( k «> k '«> k « ~ k '« - q« + q'« + q"uH(q«> q'u> q"«) x 

</3,k'u,q'u,q"u 



X 



SW(|qu — q'u — q ul)Ck'uCq'«Cq"'uCku-k'u-qu+q'u+q"u^(T ~ t'). 



Suppose, that the space of the system under study is two-dimensional. If in this space 
there is only one unstable mode with the wave vector k c and amplitude £ kc the equation 
i fTTj) is significantly simplified 
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Or 







kc 



« kc + b£ c )w - (c + d£ c + e& c ] 



dw 



kc 



Here 



A x (k c ) - -]T^(k c ,k c )s w (0) -Y,Mkc,KK)Of\k c )p^g w (k c ), 



b = u(k c , k c , k c , k c , 0) - ^ z^(k c , k c , k c )g w (k c 



(18) 



d = k c )<7 w (0) + x: ^(k c , k c , k c )0;( 1 )(k c )^°)^ w (k c 



e = ^( k o k o k c )^ w (k c ) 



The stationary solution of equation ( 1181) has the form: 



^ kc ) = iv|c + ^L + ^ c F 

= iV|c + deL + <kcl^exp 
Here iV is the normalization constant: 



2e£ 2 +d-Vd 2 -4ec 



2e§? +d+V<i 2 -4ec 



4e\/ _4 ec 



2ev / 4ec^d 2 



arctg 



\/4ec-<i 2 



ci 2 > 4ec, 
4ec > ci 2 . 



(19) 



+oo (, 

N=l/J \c + da c + e& c \-e 

— oo 

?2 i „ci 



2e£ 2 +d-Vd 2 -4:ec 



2e§ 2 +d+Vd 2 -4ec 



d 2 — 4e 



kc; 



V / |c + ^L + eekcl 4e exp 



-oo 



2ae—bd 



2eV4ec-d 2 



arctg 



\/4ec— d 2 



d 2 > 4ec, 
c?^ kc , 4ec > d 2 . 



(20) 



In the following section noise-induced effects, which arise during Turing pattern formation 
in the well-known biophysical system, will be studied. We compare the analytical results 
obtained using the above-developed approach with the results of numerical experiments. 



V. NOISE-INDUCED EFFECTS IN ONE COMPLICATED BIOPHYSICAL 

SYSTEM 



A mathematical model of the system under study is described by equations [381440 



dxi 
dt 

at 



rx 1 (l-x 1 )-^ i +D 1 V 2 x 1 , 



(21) 
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where X\. Xo are state functions, parameters r, a, b, m, g, h, f, Di, and D 2 are described in 



detail in [38l. |40(| . The investigation of local dynamics and bifurcation analysis of system (J2T 



is carried out in papers |39|,|40|. 

We introduce the dimensionless time r = rt and coordinates x' = x.<Jr/Di, and 
represent the parameters m/r and a/r in the form of: m/r = (mo/Vo)(l + /i(x',t)), 
a/r = (a /r )(l + /^(x', r)). Here mo, r , ao are spatio-temporal averages of the correspond- 
ing parameters, random homogeneous isotropic fields /j(x', r) determine the spatio-temporal 
Gaussian fluctuations of these parameters and have zero means and correlation functions of 
the form ([3]), (fTUl) . Taking into account the noise we obtain 



& = x 1 (l-x 1 



dxi 
& 
8x2 
dr 



2SL, 
''() 



^(l + / 2 (x',r)) T ^ r + J D 1 V' 2 x 1 



1 + /i(^t))x 2 -^/ 



D 2 V 2 x 2 



(22) 



A. Analytical research 

In this section the results of the analytical investigation of system (22) obtained on the 
basis of equation (fl5l) . (flBj) . (fT8|) - (1201) are given. 

For system (1220 a dispersion equation for the unstable mode averaged amplitudes of 
the form ffl5l) was obtained. Dependencies of real parts of eigenvalue A of ^££u) on the 
wavenumbers are shown in Fig. 1. 

It is evident from the dependencies in Fig. 1(a) (dotted and dash-dot lines), that the 
region of instability of the system (T22]) with ReX > increases with the intensity of external 
noise 9 2 in the supercritical region. The absolute values of the increments of unstable mode 
averaged amplitudes are greater in the presence of noise than in its absence. Therefore, 
unstable mode amplitudes increase, on the average, considerably faster in the presence of 
external noise than in its absence, which should accelerate the process of spatial pattern 
formation as compared to the deterministic description. In other words, the destruction of 
the homogeneous state and pattern formation should take place earlier in time. Besides, the 
expansion of the instability region and, consequently, increasing the number of interacting 
unstable modes should result in the changes of the form of the patterns. 

Fig. 1(b) shows the dependencies ReX(k) in the subcritical region. It is obvious, that 
instability does not appear in the absence of noise (solid line in Fig. 1(b)) or in case of its 
low intensity (dotted line in Fig. 1(b)). However, starting with certain intensity (dash-dot 
line in Fig. 1(b)) there appears a region with ReX > 0, i.e., there exists critical intensity 
of noise, that induced parametric instability of system (22). Phase transition with pattern 
formation in the presence of noise should take place in the subcritical region (i.e., earlier by 
the value of the control parameter than in the case of the deterministic description). 

In the subcritical region the state of system is homogeneous and statistically stationary 
(disorder). Unimodal probability density is known to correspond to this state. The appear- 
ance of inhomogeneous statistically stationary state (order) manifests itself in the splitting 
of probability density maximum into two symmetrical ones. 

The Fokker - Planck equation for the critical mode probability density of the form f fTSj) 
has been derived for the system fl22|) . Figures 2 and 3 show the variation of stationary 
probability density of this mode with increasing noise intensity when passing through the 
bifurcation point of a deterministic system. 
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Re^ 




FIG. 1: Real parts of eigenvalue A of unstable mode averaged amplitudes versus wavenum- 

bers k for system (22) when the noise intensity 62 changes, (a) Supercritical region (D=150). (b) 
Subcritical region (15=100). For the comparison, the figure shows curves for #i=#2=0 (solid line). 
Other parameters of the model are r =l, a =8, 5=1.434, /=0.093, /i=0.857, 6=11.905, m =0.490, 
r/i = r/2=l, #i=2.4-10~ 5 . The critical control parameter is D=135. 




FIG. 2: Steady-state probability density (|19p for critical mode values of system (|22p in the sub- 
critical region for two values of noise intensity 62- The system parameters are as in Fig. 1(a). 
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FIG. 3: Steady-state probability density (fT9l) for critical mode values of system (f22l) in the super- 
critical region for five values of noise intensity #2- The system parameters are as in Fig. 1(b). 
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Fig. 2 illustrates the probability density of critical modes in the subcritical region. If the 
noise is low the stationary probability density is close to the ^-function; the average and most 
probable values of the critical mode coincide and are equal to zero, i.e., the homogeneous 
statistically stationary state of system is the most probable one (see top Fig. 2). Deforma- 
tion of the curve of stationary probability density occurs if the noise intensity increases (see 
bottom Fig.2): the maximum is considerable reduced without being shifted; the base of the 
curve herewith is expanded. As the probability distribution obtained is not Gaussian the 
mean value becomes different from the most probable one. Thus, although the unimodal 
probability density is not split into bimodal density the mean value of the order parameter 
becomes different from zero and we should expect the inhomogeneous statistically stationary 
state (order) to arise, the probability of which is not great. It is obvious, that the occur- 
rence of a new state is accidental in the case under consideration. It can be explained in 
the following way. Formation of patterns in the subcritical region may occur on the strong 
inhomogeneities of medium. Such inhomogeneities may be produced by strong (large-scale) 
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FIG. 5: The boundary of the "order-disorder" transition on the plane of parameters D and #2- 



fluctuations, the probability of which is not great and depends on the noise parameters. 
Thus, given prolonged observation of the system's evolution and proper parameters of ex- 
ternal noise we might expect the formation of such random inhomogeneities, that will cause 
patterns to occur in the area of their location. 

Fig. 3 demonstrates the steady-state probability density of the critical mode (fT9|) of sys- 
tem fl22l) in the supercritical region at various noise intensity. Bimodal probability densities 
correspond to the existence of spatial patterns. Herewith, the most probable value and the 
expectation of the order parameter become different from zero. Fig. 3 clearly shows, that 
maximums gradually merge as the noise intensity increases, and unimodal density reappears 
at some critical value of noise intensity. At the same time ^ cmp = 0, and (£kc) 7^ 0. System 
(|22|) transforms to the state of strongly irregular behavior (disorder). Thus, the obtained 
variation of steady-state probability density of critical mode testifies, that there is a phase 
transition "disorder - order - disorder" in system (|22|) . 

The above-mentioned variation of the statistically stationary mean and most probable 
values of the critical mode corresponding to the variation of densities shown in Fig. 3 is 
illustrated in Fig. 4. 
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FIG. 6: Susceptibility x ( a ) an d second-order cumulant k 2 (b) as functions of noise intensity 9 2 - 



Fig. 5 presents the boundary of noise-induced phase transition "order-disorder" for sys- 
tem ([22]) . that is predicted using the approach developed above. It should be noted, that even 
very small fluctuations will contribute to the loss of stability of the inhomogeneous state and 
bring about a disordered state when we approaching to the deterministic transition point. 
Let us define the relative fluctuations of the order parameter (susceptibility) as x = 

(£kc) ~ (£kc) 2 / 02 and its second-order cumulant k 2 = (££ c ) j (£kc) 2 in accordance with (l2| . 

n Fig. 6 shows plots of x an d k 2 as noise intensity functions. The presence of susceptibility 
maxima (see Fig. 6(a)) clearly shows the increase of fluctuations in the vicinity of two critical 
points. The form of the curves k 2 (see Fig. 6(b)) and x obtained for system (122]) coincides 
qualitatively in the corresponding region with the similar curves obtained numerically in 

121 ] for system, in which a pure noise-induced transition is observed. This testifies, that the 
noise-induced transition "order-disorder" is of the same nature in system (122]) . 



B. Simulation 



This section presents the results of simulation of the evolution of system (122]) in two- 
dimensional space. 
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The following set of parameters m Q /r = 0.490, a /r = 8, g 2 /r Q = 2.056, / = 0.093, 
h = 0.857, b = 11.905, D 2 jD\ = 150, periodic boundary conditions and a rectangular 
domain of integration are chosen to model the evolution of system (22) in the vicinity of the 
Turing bifurcation point. Three plane waves with a low amplitude and critical wavenumber 
propagated at an angle of 60° to each other destabilized the spatial homogeneous state and 
provided hexagonal symmetry of the initial state. The size of the integration domain is 210.0 
x 186.8. 

The fluctuations of the parameters m/r and a/r are modeled as homogeneous isotropic 
Gaussian fields with the zero mean and correlation functions of the form fl3]), (II ip . where 
rg=l, k ti =100. The values of k ti inverse to the correlation time are chosen so, that all 
characteristic times of a deterministic problem are considerable greater than the correlation 
time. The choice of such values of k ti ensures, that the appropriate condition is complied 
with in theory. 

The process of spatial pattern formation in the vicinity of the Turing bifurcation point 
with increasing noise intensity 9 2 is presented in Fig. 7. In addition to the typical pattern of 
the state function x\ distribution over the surface the Fig. 7 also shows the top view giving 
a more visual representation of the structures configuration. To provide the top view the 
colour gradient from black to white visualized the variation of values of x\ from minimum to 
maximum, respectively. Black areas ( "cavities" ) correspond to empty regions of space. Fig. 
7 shows, that in the presence of noise patterns begin to form earlier than in the deterministic 
case. Besides, the greater the noise intensity, the earlier the process of formation starts. It 
should be noted, that noise destroys the symmetry of spatial patterns. 

From Fig. 7 we notice, that beginning with certain noise intensity the patterns formed 
become unstable; alternation of various random pattern configurations takes place; contours 
of certain "cavities" change, i.e., the system transforms into the state of irregular behavior 
- disorder. 

Fig. 8 demonstrates the evolution of spatial patterns arising spontaneously away from the 
Turing bifurcation point. The initial conditions correspond to the homogeneous stationary 
state of the system without initial perturbation. The boundary conditions are periodic. 
The domain of integration is square 200x200. The modeling parameters are m Q /r = 0.490, 
a /r = 8, g 2 /r = 2.056, / = 0.093, h = 0.857, b = 11.905, D 2 /D 1 = 1000. 

It is obviously from Fig. 8, that all the regularities described above are retained in the 
evolution of the system away from the bifurcation point. It must be emphasized here, that 
the "order-disorder" transition in this case takes place at greater values of noise intensity 
than in the vicinity of the bifurcation point. 

Fig. 9 gives a clear idea of the system's behavior at high level of noise. As opposed 
to pure noise-induced transition here the system's random behavior is conditioned by fast 
formation and destruction of a great number of patterns having random contours and arising 
in random places. 

In the light of the theory developed in our study such behavior is due to the fact, that 
the region of the system's instability considerable expands at high level of noise and pattern 
formation takes place owing to cooperation and competition of a very large number of 
unstable modes with different conditions of resonant interaction. 

We have also carried out simulation of system (22) evolution in the subcritical region. 
The following parameters mo/r Q — fJ,i — 0.510, a /r = fi 2 = 8, g 2 /r = 2.06, / = 0.093, 
h = 0.857, b = 11.905, D 2 /D 1 = 1000, #i=2.6-10^ 3 , fl 2 =2.56-10- 2 ; k n = k i2 = 0.5, the 
domain 200x200, and periodic boundary conditions are chosen for the modeling. 
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FIG. 7: The evolution of spatial Turing patterns at different noise intensities near the deterministic 
transition point. On the figure the moments of computer time corresponding to a given distribu- 
tion of the dynamic variable x\ of system (|22p are indicated. The bottom row of images shows 
statistically steady state in the case of order. 

Fig. 10 presents the evolution of X\ state function distribution for the instants of time t 
= 1, 927, 940, 950 and 1120. 

The pattern arising under parametric instability has a "solitonlike" shape (see Fig. 10 
t=1120). From Fig. 10 it is also clear, that the duration of the process of destroying the 
homogeneous state in the subcritical region (see Fig. 10 t=l-927) is considerably longer 
than that in the supercritical region (see Fig. 8 t=l-40 in the case of order). As predicted 
in section V.A, the initiation of inhomogeneous state in the subcritical region may be seen 
only as a result of prolonged observation of the system's evolution. 

Thus, the results of simulation of system (122!) evolution qualitatively fully confirm the 
theoretical conclusions of section V. A comparison of Fig. 5 and Fig. 7 makes it possible to 
establish a satisfactory quantitative correspondence between the theory and the numerical 
experiment in the vicinity of the Turing bifurcation point. 
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FIG. 8: The evolution of spatial Turing patterns at different noise intensities away from the 
deterministic transition point. See explanation for Fig. 7. 
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FIG. 9: The evolution of spatial Turing patterns in a strong noise. The figure shows a strongly 
irregular behavior (disorder) of the system (|22l) . The top horizontal row: #2=0. 16. The center 
horizontal row: #2=0.64. The bottom horizontal row: #2=1-44. #i=2.4040 -5 , r/i = rj2=l- 
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FIG. 10: Noise-induced parametric excitation of spatial Turing pattern in the subcritical region. 

VI. CONCLUSION 

In our study a theory is developed, that makes it possible to predict and analyze in detail 
different noise-induced effects occurring in open nonlinear distributed multicomponent multi- 
dimensional systems both in the vicinity of and away from the deterministic transition point. 
The prediction and analysis are based on the unified point of view of the concept of order 
parameters. 

Stochastic equations have been obtained for amplitudes of unstable modes (order pa- 
rameters), as well as dispersion equation for averaged amplitudes of unstable modes. The 
dependence of eigenvalues of unstable mode averaged amplitudes on wavenumbers, noise in- 
tensity, and correlation length has been found analytically. The Fokker-Planck equation for 
the order parameters of the systems under study has been obtained. Its stationary solution 
for a critical mode has been obtained in an explicit form. 

As our theory predicts, the increments of unstable modes change, the system's instability 
region is extended, the conditions of mode resonance interaction change, and noise-induced 
parametric instability occurs in external noise. The destruction of a homogeneous state and 
the pattern formation take place faster than in the deterministic case. Our theory points to 
the existence of a noise-induced " disorder - order - disorder" phase transition in systems of 
the type discussed. 

The advantages of the approach developed lies in the fact, that it applies to multicom- 
ponent multidimensional systems. The stationary solution to the Fokker-Planck equation 
for the critical order parameter is written in an explicit form. The approach applies to a 
wider class of functions Pij(x^ x 2 , Xm+ii •••> Xn')i ^(^1, x 2 , 771+1, r] s i) including those with 
a discontinuity of the second kind. Moreover, the approach suggested does not contain an 
arbitrary element connected with the discretization of system's continuous space. It can also 
be used in the case of noise with finite characteristic spatial and temporal scales. 

The numerical analysis of a specific system of the type considered qualitatively confirms 
the theoretical conclusions and provides a satisfactory quantitative correspondence between 
the theory and the numerical experiment in the vicinity of the bifurcation point. 
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Appendix A 

Functions o$„(k, k', k"), a% r „{\t, k', k", k'"), 4j'(k,k'), £^. y ,(k, k', k") introduced in 
equations (TTU|) : 

r (i) 

e )r t,i> 



^„(k,k,k") = &of\\,)o^(k')of\k'% 

4l-r» ( k > k ', k ", k "0 = £ ffgLA' w WO? (k')Of ) (k")Of ') (k" 



^k k ' k, ) = E^io; w (k)o^(ko, 



£^.„(k,k',k") = E^^'H^ofHkOofHk") 



Appendix B 

Function w(k u , k' u , k" u , k'" u , k u — k' u ) and others introduced in equations ffT4l) : 
ti - , (k tt , k u , k u , k u , k u — k u ) = crj^k^, k u , k u , k u ) — 
[°"ii^ (ku, k'm kg) + cr[^ ( k„, h s , k 7 u )] n) , „ w . 

"11 l K S5 K U1 K uj 



Ai(k s ) 

[ (7 i2 (ku, k ; u , kg) + cr 2 i ( k M , k s , k' M )] (2)/i ,// ,/// \ 

"11 l K S) K ui K uj 
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